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We discuss spatial fluctuations in the gravitational wave background arising from unresolved 
Galactic binary sources, such as close white dwarf binaries, due to the fact the galactic binary 
source distribution is anisotropic. We introduce a correlation analysis of the two data streams of 
the Laser Interferometer Space Antenna (LISA) to extract spherical harmonic coefBcients, in an 
independent manner, of the hexadecapole moment (I = 4) related to the projected two-dimensional 
density distribution of the binary source population. The proposed technique complements and 
. improves over previous suggestions in the literature to measure the gravitational wave background 

' anisotropy based on the time modulation of data as LISA orbits around the Sun. Such techniques, 

I however, are restricted only to certain combinations of spherical harmonic coefRcients of the galaxy 

^ , with no ability to separate them individually. With LISA, m = 2, 3 and 4 coefficients of the 

hexadecapole [l = 4) can be measured with signal-to-noise ratios at the level of 10 and above in a 
certain coordinate system. In addition to the hexadecapole coefficients, when combined with the 
time modulation analysis, the correlation study can also be used, in principle, to measure quadrupole 
coefficients of the binary distribution. 
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In addition to gravitational waves from massive black holes at cosmological distances, the Laser Interferometer 
^ Space Antenna (LISA; is expected to detect the galactic binary source background, such as close white dwarf 
, binaries, with gravitational waves at frequencies between 0.1 mHz to 100 mHz At frequencies above 3 mHz, LISA 
will resolve a large number of galactic binaries Q, while at the low frequency end, the unresolved population will 
] form a confusion background of gravitational waves. This background is expected to be anisotropic given the fact 
I ' that galactic binaries trace the density distribution of our galaxy such that the cumulative gravitational wave flux is 
O ] expected to be concentrated towards the highest density regimes such as the disk and the bulge of the Milky Way. 
^ As LISA orbits around the Sun, its response will be sensitive to different regions on the sky and one expects the data 
^ ' stream to modulate as a function of time. This modulation, in return, can be used to extract information related to 
. . , the density distribution and to reconstruct the anisotropy of the gravitational wave background ^lEIS- Note that the 
^ ^ ' binary background can be considered as both a source of noise or a signal; For example, by treating the background 
as a signal, we can attempt to obtain certain information related to the Galactic structure through anisotropics of 
the background ^ HI . On the other hand, if the background is treated as a source of noise, we might also discuss 
how anisotropics affect the data analysis and consider certain strategies to minimize any impact 0|. 

Here, we treat the background as a signal and consider the presence of anisotropics and their potential measurement 
using a correlation between the two data streams of LISA. For this purpose, following Ref. we introduce a spherical 
harmonic moment analysis of LISA data in terms of spatial inhomogeneities in the two-dimensional projected density 
distribution when integrated over the observer to source distance. We also emphasize the subtle use of a freedom 
related to data combinations that also corresponds to an effective rotation angle of LISA detectors in the detector 
plane. The technique we discuss here compliments, and improves, studies based on the time modulation of the data 
stream and allows one to extract individual coefficients of multipole moments related to the two-dimensional binary 
source distribution; Prior techniques were only restricted to combinations of coefficients related to the monopole 
{I = 0), quadrupole (/ = 2), and the hexadecapole (Z = 4). Considering LISA observational parameters, we show that 
certain coefficients of the hexadecapole multipole moments can be extracted with signal-to-noise ratios at the level of 
10 and above. 

We also study the complimentary nature of the two methods involving time modulation and the correlation of data 
streams and suggest that, in combination, one can also, in principle, extract information related to coefficients of the 
quadrupole distribution. The only multipole moment coefficients that will remain individually undetermined are ones 
of I = 0,m = mode (monopole) and the I = 2,m = mode related to the quadrupole. These two, however, can 
be determined, in combination, as a sum. The correlation analysis would be one of the crucial elements for LISA's 
follow-on missions aiming to detect weak gravitational wave background from the early universe . While it is 
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expected that the basic design of any potential missions in the long term is likely to be largely different from that of 
LISA, the correlation analysis with LISA data streams may present an important first step for any planned studies 
in the future that attempt to use similar procedures to both extract signals or to remove foreground noise. 

The discussion is organized as follows: in the next Section, we outline the spherical harmonic formulation of the 
anisotropics in the gravitational wave background and how the low order moments, ^ = 0, 2 and 4, are related to the 
LISA data streams. In Section 3, we discuss the correlation between data streams as a way to extract information 
related to these anisotropics beginning with a brief discussion of prior proposed techniques related to modulations in 
the data streams. In Section 4, We conclude with a summary of our main results. For background information related 
to the technique presented in this paper, we refer the reader to Ref. Q. 

II. FORMULATION 

In this section, we first discuss the measurement of gravitational wave background with LISA detectors briefly and 
then discuss how this data is related to the density distribution of the binary background, in terms of the spherical 
moments of inhomogeneities in the projected distribution on the sky. For the purpose of this discussion, we make use 
of the signal matrix involving the two data streams. 

A. Detector Response 

At the low frequency regime of the LISA band, two separate modes, {A, E), in gravitational waves can be indepen- 
dently measured such that the laser frequency noise is reduced to a level below other detector noises. Furthermore, 
individual detector noises for these two modes do not correlate . We write these two data streams di (I = A,E), 
made with detector noise nf, when the response to gravitational waves, hj, as dj = hj + nf. We define the detector 

noise spectrum Sfj as (nf*{f)n^{f )^ = 5{f — f )Sfj{f)/2. The frequency dependence on noise is almost irrelevant 

for the present discussion and, hereafter, we omit the explicit dependence for notational simplicity. Since the noise in 
two modes are uncorrelated and have same magnitude due to the symmetric data combinations, we have the following 
relation for the detector noise spectrum matrix 

^AA = ^EE^ ^AE — 0- (2-1) 

The responses hi to the low frequency gravitational waves for two modes (/ = A, E) can be essentially regarded 
as that of two 90°-interferometers rotated by 45° from each other, as shown in Fig. 1. To describe gravitational wave 
measurements, we take a reference coordinate system {X ^ ,Y£) , Z o) from the detector system, where the Zo-axis 
is normal to the detector plane. The gravitational waves behave as a Spin-2 field such that one can make linear 
combinations from the original modes [(iyi,rf_E] with a two dimensional rotation matrix i?(0) as 

('^S:;!)^'^<^*"('^)^ '''' 

These new modes [/i^ /i are equivalent to the responses of two 90°-interferomctcrs that arc obtained by 
rotation of the original modes {Hai with an angle 0/ around the Z^i-axis. Based on this transformation, one can 
easily confirm following relations for the detector noise spectrum for an arbitrary rotational angle 

SiA{<t^f)^S^E{<^f)^SiA = S'iE. S^eIM^O. (2.3) 

As we discuss later, this freedom related to the effective rotation (or data combination) is important for the measure- 
ment of anisotropics related to the galactic gravitational wave background. 

The signal matrix Sfj for the galactic binary background is defined as in the case of detector noise and is given by 
the three-dimensional spatial distribution of gravitational-wave emitting binaries, p, when integrated over the volume, 
as 

Sfj = P j d^rp{r)r-^Fij{e, 4,) , (2.4) 

where P is a normalization factor that is not important for our analysis. The weighing factor r^"^ in the integral 
is related to the fact that the gravitational wave amplitude scales as r~^, where r is the radial distance from the 
detector. Here, we assume that the frequency distribution and the spatial distribution of galactic gravitational wave 
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emission rate are independent. The functions F/j represent the dependence of the detector's response on the source 
direction, and are defined as 



1 + cos^ 



(2.5) 



with {pi,p2) = (cos^ 2^, sin^ 20) for Faa, iPi,P2) = (sin^ 2(?!), cos^ 2(/)) fo r F ee, and (pi,_P2) = 
(sin 2(/) cos 2(^, — sin 2(/) cos 20) for Fae- Using a spherical coordinate we can express Eq. (E3I) as 



sfj^ / dnB{e,ci))Fij{e,<t>), 



(2.6) 

where the angular dependence of the gravitational wave background luminosity B(6, 0) is defined through the radial 
integral as 



P j drp{r,e, 



(2.7) 



Note the dependence on factor in Ea. l2.4l is removed due to the r^dr factor involved in the d^r integral. Thus, the 
projected two dimensional structure of the density distribution, as seen by gravitational waves, is simply proportional 
to the column density of the binary population; this is also proportional to the radial distance averaged density profile 
and, roughly speaking, the function B(9, (j)) can be regarded as the optical luminosity distribution on the sky. As 
before, one can relate the signal matrix under a rotation such that Sfj{<j>f) = R{2<j) f) Sfj R{—2(f> f) for the new data 
combination of [dA{<j>f),dEi4>f)]- 



B. Spherical Harmonic Expansion 

As shown in Fig. 1, the noise matrix Sfj is determined by the orientation of the coordinate system Kd', {Zd, Yd, Zd) 
which changes with the rotation of LISA. Given this, we consider the spherical harmonic expansion under a fixed 
coordinate system Kq; {Xq, Yq, Zq) since such a fixed system can aid in analyzing the gravitational wave background 
In this system, we set the {Xq, lo)-plaiie on the ecliptic and take X^-axis in the direction of the autumn equinox 
from the Sun. This is a reasonable choice considering the LISA configuration (l|, [7| . The two systems Kjj and Kq are 
related by the Euler angles (a,/3,7). Here {a, (3) is the direction of the Z^-axis in the fixed Kq system, while 7 is 
essentially degenerate with the freedom related to (/>/-rotation in Ea. H2.2|l around the Zu-axis. Therefore, we use the 
angle 7 hereafter as an independent variable that can be varied. 

The angular dependence of the gravitational wave intensity -6(0,0) can be decomposed with spherical harmonic 
moments Y/„i(0, 0) (with the angular variables in the Kq system) as 

B(0,0) = ^B,„,yz„(0,0). (2.8) 

l,m 

Here, harmonic coefficients Bi^ are given as 

Bira = {lm\B) , (2.9) 

where we have used the traditional notation for the inner product {a\b) = j dna*{9,(f>)b{9,(j>), and an abbreviation 
I^to) = \Yim{0, 0)). Our goal in this subsection is to write down the matrix Sfj{a, /3, 7), in the moving Kd frame that 
is characterized by the Euler angles (a,/?, 7), with coefficients Bim related to the background density distribution. 
First, the matrix Sfj{a,f3,j) is formally expressed as 

SPj{a,(3,j) = (F,j(0,0)|C/(a,/3,7)"i 15(0,0)) , (2.10) 

where C/(q;,/3, 7) is the rotation operator related to Euler angles a,f3 and 7. After some tedious but straight forward 
algebra (see, the Appendix A of Ref. for details), we write 



\-ma Izl 



SEe = ^Soo + 7? E KnMe-'-B,,,,] + ^ 5: i?o^„Re[e™"S4™] " M E DtMB.n.e^'''^-^'] 



^AE = ^J2D'^r.M[B,me^""-'-^']- (2.11) 
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As written, the two data streams from LISA detectors are only sensitive to the background Bi„i with I = 0,2 
and 4 harmonics. This is simply due to the angular dependence of the Fij{9,4>) term such that {Fjj\lm) is only 
non-zero when ^ = 0, 2 and 4. Unlike electromagnetic radiation observations, such as say that of photons in the 
cosmic microwave background, where, in principle, all multipole moments of the inhomogeneous radiation field is 
measurable, with gravitational waves, one is only restricted to these low order moments of the underlying field due 
to the poor directional characterization of the gravitational wave detectors. One can potentially extract additional 
multipole moments of the gravitational wave background with more data combinations involving, say, an additional 
configuration such as another set of 3 spacecrafts . Since the LISA configuration is fixed to three spacecrafts, 
we do not discuss such possibilities further. Note that our expressions in Ea. (|2.11() are valid at the long wavelength 
limit of LISA relevant for the Galactic confusion background. At higher frequency regime, we expect the dependence 
of these terms in terms of harmonic coefficients to change 0, 0| • 

The off-diagonal element of the noise matrix S^^ is non-zero due to the I = 4 coefficients of the binary background. 
The above equations related to Sfj{a, (3,^) are essentially same as eqs.(A13)-(Al5) in Ref. 7]. Here, in comparison, 
we have rewritten the spin-weight harmonics in that paper with matrix elements D[j^{(3) — {^ls\U{0, (3,0)~^\lni) with 

-Mf3,a)^l^^ = DUP)e^^- (2.12) 

and elements D^^dS) given by the Jacobi polynomials p^(^+™''' When s>m, their explicit forms are 

D\^{(5) = (-1)' 

while for s < m, Di^{f3) = (-l)'^-™7^5„,(/3) 

In the case of observations with LISA, the angle /3 is fixed to the value of 7r/3. This leads to the restriction that, 
in principle, one cannot separate information related to Bqq and 820- The angle a is given as a = 27rt with orbital 
time t in units of year (measured from the Autumn equinox). For the Galactic binary distribution p(r) we use the 
triaxial model given in |0 and used in f?\. The relevant coefficients Bim and -D^,„j are presented in Tables I and II. 

While we primarily discuss the background from our galaxy, it is worth considering the possibility that anisotropics 
can exist due to gravitational wave emission from binaries in other galaxies of the Local group, such as the Andromeda 
galaxy (M31) or the Large Magellanic Cloud (LMC). Their contribution to the anisotropy B{6,(j)) is roughly pro- 
portional to (mass)/(distance)^. The characteristic distance is lOkpc for the Milky Way galaxy, 50kpc for LMC and 
TOOkpc for M31. As for the mass, M31 is comparable to Milky Way, but LMC is ~ 50~^ times smaller. Given these 
parameters, we find that the contribution to anisotropics of the gravitational wave background from other member 
of the local group is negligible. This is also apparent from the fact that the optical luminosity of these and other 
galaxies are much lower than the Milky Way. 

III. EXTRACTION OF ANISOTROPY PARAMETERS 

Having established the relation between components of the LISA data stream matrix and coefficients of the binary 
background anisotropy, we now discuss the individual measurement of these anisotropy coefficients based on various 
combinations of terms in the data matrix. First, we briefly outline the previous method based on the time modulation 
of the data stream and consider the relation between anisotropy coefficients there and the ones we discuss here. 

A. Time Modulation of LISA Data 

Previous studies on the measurement of anisotropics Bim are based on the time modulation of the spectrum 5'_f^ 
0, 0, 01 (see also 0). These techniques do not use (i) other modes S^^ or and (ii) the freedom to adjust 
one of the Euler angles, 7, a priori. In time modulation studies, the angle 7 simply follows the motion of LISA and 
given as 7 = —2Trt + 70 where 70 is some constant angle (here, we put 70 = for simplicity). We briefly follow the 
underlying idea behind time-modulation method and highlight certain problems for the signal extraction related to 
The function S"!^ can be expanded as Sf^{t) = S^^{a = 2nt,f3 = 7r/3,7 = -2Trt) = EL-8 ^'"^2'^™*. Since 
C* = C-n, as one is observing a real realization on the sky, we do not discuss C„ with n < 0. The explicit forms of 
the coefficients C„ are easily derived from Eas. H2.11|l 

Co = -|— -Soo + ^^^-000^20 + -^40^00 + -^]J^Dj^Re[Bi4], (3.1) 



{i~sy.{i + sy. 

{I - my.{l + m)l 



Us|j fsinfj p/^r"-™)(-cos/3) 



(2.13) 
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FIG. 1: Definition of the detector coordinate {Xd, Yd, Zd)- The Zo-sxis is normal to the detector plane. A and E modes can 
be regarded as two L-shaped detectors rotated by 45°. 



5.5. X 



Galactic center 




Autumn equinox 



Ecliptic plane 



FIG. 2: Configuration of LISA and the definition of the fixed ecliptic coordinate {Xo,Yo, Zo). The ^o-axis is normal to the 
ecliptic plane and Xo-axis is oriented to the autumn equinox from the Sun. 
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(3.3) 
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(3.4) 
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Assuming that the detector noise S^j^ is significantly smaller than the background S*^^ in a frequency bandwidth 
A/, the signal-to-noise ratio for each coefficients (n 7^ 0) is given as 



\Cn\ 



(3.6) 
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where T is the observational period. Statistics for the measurement of the time dependent signal C„ (n ^ 0) would 
be largely different from that of the constant part Cq. The coefficients C„ are complex numbers but their phases 
depend on the choice of the orientation of the Xg-axis that is somewhat arbitrary in the ecliptic plane. Therefore, 
we summed up the signal to noise ratios of their real and imaginary parts in Eq. 1)3. 6|) . The numerical prefactor of 
•\/2 comes from the fact that we assume the noise contributions to the real and imaginary parts of these coefficients 
are related to the total noise, Nt such that Nr ^ Nj = Nt/\/2 without correlation between them. For our model 
parameters related to the binary distribution, we have 

(3.7) 

with T = lO^sec and A/ = 10 ''Hz. Thus, one can obtain only Ci and C2 with sufficient statistical significance. 
The above signal-to-noise ratios are consistent with previous analyses when considering the difference between models 
for the Galaxy density distribution between these calculations and the present one Since Ci and C2 depend on 
certain combinations of Bim, one cannot separate out the multipole moment coefficients related to the gravitational 
wave background by simply measuring the time modulation components in the Fourier space. Moreover, in terms 
where certain Bim coefficients are simply proportional to C„ coefficients, especially when n > 5, the dependence is 
significantly reduced through small coupling coefficients. For example, the information related to -B4±3 appears in 
C7 with the factor —Df_^y^Tr/70/3 — —0.0044. Once can consider combinations of C„ to extract a certain Bim, but, 
again for example for the case with i?4±3, the dependence in C3 is through DQ^^/^^/105 ~ —0.0081 and in Ci through 
— Z3|3-\/7r/70/3 = —0.036. The former factor (in C7) is smaller than the latter (note -D|_3 ^ D^q, Df^^; Table II), but 
information related to in the latter remains to be degenerate with other Bim coefficients such that one cannot 
easily find the necessary combinations. The situation is also similar for i?4±2 or i34±4. The factor £'|_4^7r/70/3 in 
Cg is extremely small, and it is impossible to determine i?44 form the coefficient Cg. 

From these examples, we learn two things: We should try to improve the signal-to-noise ratio for Bim coefficients by 
taking a suitable combination with respect to the matrix element Df^, and handle the degeneracy of Bim coefficients 
with the use of other modes. 

In the next subsection, we study a method for signal analysis where we pay special attention to these two points. 



B. Correlation Analysis 

We use the correlation analysis / SAE{t)dt to measure the parameter B4,n- By setting 7 — — 7r/2m't-|-7r/8 (m' > 0), 
for each m' from to 4, a priori — since we have this freedom to specify this angle before hand — one obtains, under 
the scenario that a — 2'k t and after integrating the signal SAsit) over an observational period T (an integer in units 
of one year). 



2V7r „4 



3V70 



Dt„^,Rc[Bi^.] , (3.8) 



while the imaginary component of the same quantity can be obtained by setting 7 = — 7r/2m't such that (S'^b)™ oc 
Im[i34„i]. Thus, the substitution for 7 has the advantage that one can now individually get real and imaginary parts 
of B4m, while at the same time, also making use of the substantially large coupling coefficients of Df^^ (m > 0) that 
are associated with harmonic coefficients. 

When we neglect the detector noise and the frequency dependence of anisotropics, which would be tolerable for the 
present study, the signal-to- noise ratio of the correlation analysis, for an individual I — 4 and m coefficient, is 




(S/N)4™ ^\ dt df '' ^^fj '^^ , (3.9) 



where we have now combined the real and imaginary parts [{SaeY"]'^ = [('^'ab)^]^^ [("^ab)™]^' have numerically 
calculated these signal-to-noise ratios for each m value using the same density distribution model for the galactic 
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l,m 


0,0 


2,0 


2,1 


2,2 


4,0 


4,1 


4,2 


4,3 


4,4 


Blm/Boo 


1 


-0.243 


-0.292+0.077i 


-0.647-0.038j 


0.052 


-0.076+0.076i 


0.289-0. 012i 


-0.329-0.073i 


0.504-h0.061i 



TABLE I: Spherical harmonic coefficients of the Galactic binary background. Note that Bj-m ~ (~1)'"A* 















A^ 






0.459 


0.530 


-0.125 


-0.530 


0.459 


^0,-4 




^0,-2 


At,-i 


^0,0 


^0,1 


^0,2 


^0,3 


^0,4 


0.294 


0.480 


0.222 


-0.303 


-0.289 


0.303 


0.222 


-0.480 


0.294 










Dio 




^?l2 


AS 


i3t.4 


0.004 


0.019 


0.062 


0.152 


0.294 


0.456 


0.558 


0.516 


0.316 



TABLE 11: Matrix elements D[^{I3) for (3 = 7r/3. 



binaries as before and we obtain following values for {S/N)im for the measurement of B^m coefficients, separately, as 

(5/iV)4o = 1.1, (5/iV)4i = 3.7, (5/iV)42 = 12, (5/iV)43 - 13, (5/A^)44 = 12, (3.10) 

when A/ = 10~'^Hz and T — lO^sec. Thus, as estimated, we expect one can estimate the coefficients i?42, i?43 and 
^44 accurately, while upper limits can be established for B^q and i?4i. Then, using the previous time modulation 
method, one can further estimate B21 form Ci coefficient and B22 from C2 coefficient. Since coefficients Ci and C2 
are dominated by i3i2 and ^22 respectively, this determination can be carried out with sufficient signal-to-noise ratio. 

So far, we have studied the correlation analysis for estimating i?4m individually for all m values and the combination 
of time modulation related to S*^^ and the correlation analysis to extract certain the coefficients related to the 
quadrupole. One can also use the time modulation method alone, but with combinations of S'^yj, S^^^ and S^^. An 
example is the modulation of 5'^^ — S^^. Again, the freedom related to the angle 7 is important. For example by 
setting 7 = 27rt(— m/4 -I- 2) or 27rt(— m/4 -I- 2) -I- 7r/8, we can estimate the real and imaginary parts of -64™ from the 
coefficient Cg related to the combination 5*^^ — S^^. By making numerical estimates, we find that the signal-to-noise 
ratio for such a measurement is similar to the one estimated above with the correlation analysis. 

While the combination with time modulation for alone do lead to certain coefficients of i?4m, we find that the 
correlation analysis improves the correlation method by providing a mechanism to extract all coefficients. Moreover, 
the correlation method compliments techniques related to the time modulation and we find that, in combination, the 
methods can be used to extract most coefficients related to the / = 2 and 4 moments of the galactic distribution of 
binaries. The only coefficients that are not independently measured is Bqo and B20, though the two can be measured 
in combination as a sum. As mentioned earlier, signal analysis for the Cq mode is different from other C„ modes with 
n ^ 0. For this purpose, data streams other than A and E modes, e.g. Sagnac, would be important to calibrate the 
detector noise (T3 |. 

IV. SUMMARY 

In this paper, we have discussed spatial fluctuations in the gravitational wave background arising from unresolved 
Galactic binary sources, such as close white dwarf binaries, due to the fact the this galactic source distribution, 
when projected on the sky, is anisotropic. We introduce a spherical harmonic analysis of this anisotropy based on 
a correlation study of the two data streams of LISA and propose the measurement of individual spherical harmonic 
coefficients related to the hexadecapole moment (Z = 4), in addition to coefficients of the quadrupole (/ = 2). 

The proposed technique complements and improves over previous suggestions in the literature to measure the 
gravitational wave background anisotropy based on the time modulation of the single data stream [e.g. S^j^{t)) as 
LISA orbits around the Sun and restricted only to combinations of spherical harmonic coefficients with no ability to 
separate them. The correlation analysis, by making use of a single freedom related to the rotation angle of detectors 
in the detector plane or similarly how data are combined, provides a method to independently extract multipole 
coefficients of the hexadecapole moment of the background. The information related to the quadrupole comes from 
the correlation analysis in combination with the time modulation. The only independently undetermined coefficients 
turn out to be that of the monopole and the to = mode of the quadrupole, though the two can be determined as a 
sum. 
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With LISA, we have shown that certain coefficients of the hexadecapole can be measured with signal-to-noise ratios 
at the level of 10 and above. While LISA will present a first step towards a correlation analysis of the gravitational 
wave background, future missions are likely to exploit this in detail both as a source of signal, to extract spatial 
information, and noise, to decrease noise to detect almost isotropic backgrounds such as the one related to primordial 
gravitational waves. This would certainly be the case at low frequencies as the foreground is likely to be the dominant 
source of anisotropy. At few mHz frequencies, the fluctuation amplitude related to the extragalactic component is 
expected to be at the level of a few percent, but restricted to arcminute angular scales where galaxy density field is 
observed to vary. At /iHz frequencies, where the merging massive black holes are likely to be the domiant foreground 
source, the anisotropy will be dominated by the shot-noise related to the finite number density of such binaries on the 
sky. On the other hand, the anisotropy of the primordial gravitational wave background is expected to have a pattern 
close to that of the cosmic microwave background, with a fluctuation level of <5 ~ 10~^, at large angular scales. A 
confirmation of their similarity in anisotropy can be used as a proof on the primordial nature of the background. For a 
direct detection of spatial fluctuations in the gravitational wave background, in the optimistic case where the detector 
noise is dominated by the monopole mode of the background, one is forced to gravitational wave frequencies at the 
level of / ~ A/ - 10''(SN/10)2(,5/10-5)-2(r/3yr)Hz. Thus, for a detection of the anisotropy at a level S - IQ-^ 
over a realistic observational period with significant signal to noise ratio, the relevant frequencies are at / > lOkHz. 
Interestingly, at this high end of the gravitational wave spectrum, there are no viable astrophysical sources and 
suggests that, in future, one may in fact be able to extract anisotropy associated with the primordial gravitational 
wave background. 
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